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Fe and N self-diffusion in non-magnetic Fe:N has been studied using neutron reflectivity. The isotope labelled 
multilayers, [Fe:N/^^Fe:N]io and [Fe:N/Fe^^:N]io were prepared using magnetron sputtering. It was remark- 
able to observe that N diffusion was slower compared to Fe while the atomic size of Fe is larger compared to 
N. An attempt has been made to understand the diffusion of Fe and N in non-magnetic Fe:N. 
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Iron nitrides (Fe:N) show a variety of structures and 
magnetic properties with a variation in the nitrogen con- 
tent. With an increasing atomic percentage {at.% N), 
their major phases are: Fei6N2, Fe4N, FesN, Fe2N, 
FeN and Fe3N4. For <25 at.% N, the Fe:N phases 
are magnetiCfi A lot of attention has been driven to 
a" — Fei6N2 (~11 at.% N) due to the presence of the 
so-called giant magnetic moment in this compound. ^"'^ 
Around 20 at.% N, 7' — Fe4N phase is formed which has a 
well-defined magnetic properties and crystal structurei^ 
Very recently the 7' phase has received a lot of inter- 
est due to its chemical inertness and mechanically hard 
surfaces making it a suitable alternative to pure Fe in 
magnetic deviceSf4"= Between 25-33 at.% N the Fe:N are 
known as e — Fe^jN (2 < a; < 3), and as Nat.% increases 
from 25% to 33%, the phase changes from ferromag- 
netic FeaN to paramagnetic Fe2N at room temperature. 
Earlier it was not possible to produce the Fe:N phases 
containing more than 33 at.% N but by using reactive 
sputtering^"^^ and pulsed laser deposition techniques, 
the Fe:N phases with >33 of.% N were produced. Around 
50of.% N, the FeN phases have cubic ZnS and/or NaCl- 
type structures which are known as 7" and/or 7"'. The 
Fe3N4 phase with even more than 50 at. % N was pre- 
dicted by Ching et alJ^, but has not been evidenced 
experimentally. 

Recently, the e — Fe2N and y /y phases have been 
used as a precursor to prepare 7' — Fe4N phase for its 
use in the spintronic devicesi^^— In the present work, 
we have prepared single phase e — Fe2N and 7"' — FeN 
compounds and studied the self-diffusion of Fe and N. 
A proper understanding of the stability and nitride for- 
mation requires the knowledge of both Fe and N self- 
diffusion at atomic length scales. However, there are no 
studies on measurements of both Fe and N self-diffusion 
in non-magnetic Fe:N. Conventional techniques to mea- 
sure self-diffusion (e.g. secondary ion mass spectroscopy, 
radioactive tracer etc.) have depth resolutions of several 
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nm. Therefore, to measure diffusion at nanometer length 
scale, a technique with depth resolution in the sub-nm 
regime is necessaryi^ Here, the method of choice is neu- 
tron reflectometry (NR) , which in addition is sensitive to 
isotopic contrast: neutron scattering length for natural 
Fe and ^^Fe are 9.45 and 2.3 fm, and for natural N and 
^^N 9.36 and 6.6 fm, respectively. 

The samples were prepared at room temperature 
by DC magnetron sputtering using either a mixture 
of Ar and N2 (samples A), or with pure N2 (samples 
B) as sputtering gas. During the deposition the total 
gas flow was kept constant at 10 standard cm'^/min 
(seem) and the sputtering power at 50 W. The actual 
thicknesses (obtained using NR) and gas flow parame- 
ters are; Samples (A): [Fe:N(7.5 nm)/^^Fe:N(4.5 nm)]io, 
[Fe:N(6.4nm)/Fe:^^N(3.2nm)]io multilayers deposited 
using sputter gas (N2 -|- Ar) each at 5 seem. And 
samples (B): [Fe:N(10.2nm)/5^Fe:N(5.2nm)]io, 
[Fe:N(10.6nm)/Fc:i^N(5.2nm)]io multilayers using 
10 seem N2 as sputter gas. The multilayer with ^''Fe was 
deposited by alternatively sputtering Fe/^^Fe enriched 
targets and the ^^N multilayer was prepared using the 
same Fe target but switching between N/^^N gases. A 
residual gas analyzer was installed in the sputtering 
chamber to monitor the isotope abundance of nitrogen. 
The samples were characterized using x-ray diffraction 
(XRD) using Cu-Ka x-rays while thermal stability was 
examined using differential scanning calorimetry (DSC). 
The local environment of ^^Fe atoms was probed using 
conversion electron Mossbauer spectroscopy (CEMS). 
The NR measurements were performed at the reflec- 
tometers AMOR at SINQ/PS^i^ and D17 at ILlJ^, 
both in the time-of-flight mode. 

The CEMS and XRD patterns of samples A and B 
are shown in fig. [T] (a)-(d). The XRD pattern of sam- 
ple A (fig. [TJs) , shows peaks corresponding toe — Fe^N 
with hexagonal closed pack structure. The CEMS spec- 
trum of this sample (fig. [1^) shows an asymmetric dou- 
blet which was fitted using two doublets corresponding to 
Fe-III and Fc-II sites with isomer shifts of 0.43±0.002 and 
0.37±0.002mm/s; quadrupole splitting of 0.24(±0.002) 
and 0.50±0.004mm/s and the relative area ratio of 69:31, 
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FIG. 1. (Color online) CEMS pattern of sample A (a) and 
sample B (c). XRD pattern of sample A (b) and sample B 
(d). 
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FIG. 2. (Color online) Neutron reflectivity pattern (reflec- 
tivity R, multiplied by for better visuality of Bragg peak) 
of sample B measured at 463 K for different annealing time 
taken with an interval of 30min. The time dependence of Fe 
(a) and N diffusivity (b) can be seen. 

respectively. The fitted parameters matches well with the 
reported values and using the relative area ratio, we ob- 
tain the value of a; = 2.23 following the procedure given 
in ref.^ 

The XRD pattern of the sample B (fig. [TJl) shows all 
the peaks corresponding to 7"' — FeN phase with lat- 
tice parameter a = 0.454 ± 0.001 nm. The CEMS pat- 
tern of this sample (fig. [TJ;) shows an asymmetric dou- 
blet and was fitted using two singlets with isomer shifts 
of 0.05 ± 0.002 and 0.6 ± 0.003 mm/s. The XRD and 
CEMS parameters match well with the values obtained 
by Jouanny et alM- for 7"' — FeN having ZnS-type struc- 
ture with Ri50at.% N. The average grain size calculated 
using the Scherrer formula for the most intense peak in 
sample A and B arc 16 nm and 5nm, respectively. 

The thermal stability of the samples was studied using 
XRD and DSC. In the temperature range of 373-573 K 
there was no appreciable change in the XRD patterns and 
DSC measurements showed a strong exothermic peak at 
650 K, indicating out-diffusion of nitrogen as observed 
in an earlier study. Therefore 523 K was chosen as the 
maximum annealing temperature for the diffusion mea- 
surements. The ^''Fe and ^^N periodicity gives rise to 
Bragg peaks in the NR pattern, which become less sharp 
as the multilayers are annealed. Some representative NR 
patterns for sample B at 463 K at different annealing time 
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FIG. 3. (Color online) Time dependence of Fe and N instan- 
taneous diffusivity for sample B at 463 K (a). Arrhenius be- 
havior of Fe and N diffusivity for sample B (b) and for sample 
A (c). The correlation between InDo and E (d) for crystalline 
(□) and amorphous (Q) alloys^. The star (-^) represents 
data obtained in the present work and (-I-) corresponds to the 
values obtain in crystalline thin film multilayers.— 

TABLE I. Parameters for Fe;N samples. 



Sample 



A 



B 



Composition 
Epe (eV) 
EN(eV) 
InDop^(mVs) 
InDojv (m^/s) 



e - Feg.gsN 

0.25±0.03 

0.46±0.08 

-44±1 

-41±2 



7"' - FeN 

1.0±0.2 

1.5±0.3 

-25±10 
-16±10 



(with a step of 30 min) are shown in fig. [2] (a) and (b) for 
Fe and N contrast. The pattern was multiplied with 
to remove the decay due to Fresnel's reflectivity. The 
decay of the Bragg peak intensity can be used to calcu- 
late instantaneous diffusivity (Di) using the expressioni^^ 
ln[/(i)//(0)] = -Sn'^n'^D.t/cf where /(O) is the intensity 
of the n"* order Bragg peak at time t = (before anneal- 
ing) and d is the bilayer thickness. As the structure tends 
to relax, the time averaged diffusivity is defined byj^ 
D — J Di{t') dt' . Assuming an exponential law for re- 
laxation, Di{t) = Const. exp(— t/r) + D, where D is dif- 
fusivity in the relaxed state and r is the relaxation time. 
Using this relation the diffusivities in the structurally re- 
laxed state were obtained and are shown in fig. |3](a). 
The temperature dependence of diffusivity obtained in 
the structurally relaxed state follows an Arrhenius-type 
behavior as shown in fig.[3](b) for sample B. Similarly, for 
sample A, the Arrhenius behavior for Fe and N diffusiv- 
ity is shown in fig. [21(c). The straight line fit to the data 
obtained using: D ~ Dq exp{—E/kBT), yields the pre- 
exponential factor (Dq) and activation energy (E) which 
are given in table |T] with /cb being Boltzman's constant. 
It is known that a correlation exists between ItiDq and 
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E, which is followed in all class of materials and indicates 
the involved diffusion mechanism. While including our 
values of E and InDo in this correlation (fig. [3jl), we 
find that our values follow this correlation well for dif- 
fusion in amorphous alloys. With a straight line fit of 
the values obtained for amorphous alloys with our data, 
we can calculate the entropy (AS*) for diffusion following 
the procedure given in ref.^** The calculated values of AS" 
for Fe and N diffusion are 5 fee and 19 fcs for sample A 
and 9 fce and 28 fce for sample B, respectively. In crys- 
talline systems the value of AS* is typically (3-5) fee corre- 
sponding to interstitialcy or monovacancy mechanism;^^ 
The relatively large values of AS indicate that the dif- 
fusion mechanism in the present case is more similar to 
amorphous alloys where collective type diffusion occurs 
involving a group of atoms. As we see from the XRD 
measurements, the grain sizes are rather small (16 nm 
and 5 nm in sample A and B). If the diffusion takes place 
from the grain boundaries which are not so well-ordered, 
the diffusion mechanism will lead to a situation similar 
to amorphous alloys. In an earlier study, we measured 
the Fe self-diffusion in a nearly equiatomic amorphous 
and nonmagnetic FeN alloy^^ and found EFc=1.3±0.2eV 
which is close to EFc=l±0.2eV for sample B. Therefore 
nearly similar values of E in amorphous and crystalline 
FeN thin films indicate a similar type of diffusion mech- 
anism which is independent of the structure. 

The comparison of diffusivity in samples A and B re- 
veals that the values of E are greater in sample B. This 
indicates a slower diffusion of Fe or N in sample B as 
compared to sample A. Looking at the enthalpies of 
formation (A7J°), the values for FeN are the lower as 
compared to Fe2N or other Fe:N phasesi^ This means 
that FeN is expected to be more stable as compared to 
Fe2N and therefore E for Fe and N self-diffusion should 
be high in FeN as compared to Fe2N. Though the differ- 
ences in the obtained values of diffusivity in both sam- 
ples can be understood in terms of energetics of iron ni- 
trides, it is the difference between the Fe and N self- 
diffusion which is counter intuitive as the atomic size of 
Fe (rpc = 0.1274 nm) is larger than N (rpc/^'N ~ 1-6). 
Therefore, it is expected that N should be diffusing faster 
than Fe as was observed in a N-poor magnetic Fe:Ni^ In 
the literature we wee that atomic size dependency (i.e. 
smaller atoms diffuses faster) is observed metal-metal 
and metal-metalloid amorphous alloysS2j2^. However, in 
case of phosphorus self-diffusion in Fe4oNi4oPi4B6 and 
in Pd43Cu27NiioP2o metallic glass the diffusion coeffi- 
cient of P was found to be smaller compared to Fe and 
Cu i^^i^^ It was suggested that the local chemical interac- 
tions around P (strong covalent bonds) are more impor- 
tant for the diffusion than the atomic size dependence of 
constituent elements in the alloy. In the periodic table 
N and P are in the same group (VB) and the anomaly 
observed for P self-diffusion can be extended to N self- 
diffusion in nonmagnetic Fe:N phases. However, as the 
temperature increases we observed that the difference be- 
tween Fe and N diffusivity becomes smaller. 



In conclusion, we measured Fe and N self-diffusion 
in non-magnetic e — Fe2.23N and 7"' — FeN compounds. 
The activation energy of Fe and N is higher in 7"' — FeN 
than in e — Fe2.23N. However, the diffusivity of N is 
smaller compared to Fe in both samples. The differences 
in Fe and N diffusivity can not be understood in terms 
of atomic seize of Fe and N. 
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